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Abstract 



This second paper, on spin-0 amplitudes, is based on the underlying 
results and methods outlined in the preceding Paper I, which describes 
^ , the complex approach to quantum amplitudes in black-hole evaporation. 

' The main result of the present paper is a computation of the quantum 

^S) I amplitude for a given slightly anisotropic configuration of a scalar field </!> 

, on a space- like hypersurface Ef at a very late time T, given also (for 

' simplicity) that the initial data for gravity and the scalar field at an 

initial surface E/ are taken to be spherically symmetric. In particular, 
. this applies to perturbations of spherically-symmetric collapse to a black 

hole, starting from a difi'use, nearly-stationary configuration, where the 
^ , bosonic part of the Lagrangian consists of Einstein gravity and a massless, 

O^' minimally-coupled real scalar field <j!> . As described in Paper I, Feynman's 

+ie approach is taken; here, this involves a rotation into the complex: 
{3/), T \T\ exp{—i9), with < 9 < n/2 . A complex solution of the classical 

boundary- value problem is expected to exist, provided > 0; although for 
9 = (Lorentzian time-separation) , the classical boundary- value problem 
■ is badly posed. Once the amplitude is found for 9 > , one can take the 

' limit 6* 0+ to find the Lorentzian amplitude. The paper also includes 

a discussion of adiabatic solutions of the scalar wave equation, and of the 
dimensionality of certain quantities which occur frequently in the project 
as a whole, for which these two initial papers establish the underlying 
results. 



1 Introduction 

In this treatment of spin-0 quantum amplitudes in black-hole evaporation, a 
concrete calculation is given, based on the general procedure outlined in the 
previous Paper I on the complex approach to such amphtudes [1], and in [2]; 
see also [3]. Consider asymptotically-flat boundary data {hij , and {hij ,4')f, 
given on two space-like boundary hypersurfaccs S/ and E p , which are separated 
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by a large Lorentzian proper-time interval T, as measured near spatial infinity. 
Here, {hij)i^F = {gij)i^F = 1,2,3) gives the intrinsic Riemannian 3-metric 
on the surfaces E/ and T,f , where g/^i, = 0,1, 2, 3) gives the space-time 4- 
metric. As in Paper I, we assume, for the sake of dcfinitcness, that the bosonic 
part of the Lagrangian contains Einstein gravity with a minimally-coupled real 
massless scalar field <p , as posed above in the boundary data. One would like 
to calculate the quantum amplitude corresponding to these data. It is certainly 
impossible to compute this amplitude directly via a semi-classical expansion, 
that is, as approximately exp(i5ciass), where 5ciass is the Lorentzian action of 
a classical solution of the coupled Einstein/scalar field equations, subject to 
the boundary data above; typically there is no such classical solution, since the 
boundary- value problem for hyperbolic equations is badly posed [1,4,5]. Instead, 
following Feynman's +ie prescription [6] , we rotate the asymptotic time- interval 
T into the complex: T \T\ exp{—i6), for < < ■7r/2 . If the resulting com- 
plex boundary- value problem, up to gauge, is strongly elliptic [7], one would 
expect existence and uniqueness properties as for real elliptic equations. In 
particular, this should give semi-classically a quantum amplitude proportional 
to exp(i5ciass) = exp(— /class), where I denotes the Euclidean action, provided 
that the Lagrangian is invariant under local supcrsymmetry - it appears that 
meaningful quantum amplitudes can only be defined in that case [5,8,9]. Fol- 
lowing Sec. 1. 2, one would expect such a semi-classical amplitude to have finite 
loop corrections in a locally-supcrsymmetric field theory of supergravity coupled 
to supermatter, where the loop corrections will only contribute significantly for 
boundary data involving Planckian energies. For a pure supergravity theory, 
the semi-classical amplitude is expected to be exact [5,8,9]. Finally, to obtain 
the amplitude for a Lorentzian time-separation T , one takes the limit of the 
amplitude as ^ ^ 0+ . 

The concrete calculation in this paper concerns the quantum amplitude cor- 
responding to emission of scalar radiation in the black-hole collapse problem, 
as measured by non-trivial non-spherical perturbations of (j) on the final surface 
'Sp , which is chosen to be at an extremely late time T , so as to intersect all the 
outgoing radiation. To simplify matters, for the purpose of the present calcula- 
tion, we assume that there are no gravitons on Sj? , that is, that there are no 
non-spherical perturbations present in the final gravitational data hijp . Equiv- 
alently, we assume here that hijp describes an exactly spherically-symmetric 
spatial gravitational field. Later on, we shall compute the corresponding am- 
plitudes for the opposite case, for which the only non-trivial perturbations in 
the final data are gravitational, residing in the 3-metric hijp , whereas 0^ is 
taken to be exactly spherically symmetric; this gives the amplitudes for generic 
spin-2 perturbations. Further to simplify matters, we shall assume, in this pa- 
per for spin-0 and in subsequent work for spin-2, that the initial data {hij , 4>)i , 
representing the gravitational 3-metric and scalar field at an early time before 
collapse, are exactly spherically symmetric. Of course, once one has the above 
amplitudes for weak spin-0 and spin-2 fluctuations separately, one can combine 
them. 

In Sec. 2 of the present paper, we consider the {t, r)-dependent decoupled lin- 
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ear equations for the functions Ri„i appearing in the decomposition of the scalar 
wave equation into angular harmonics. Here, £ = 0, 1, 2, ... , with m subject 
to — i < m < £ , are the usual angular-momentum quantum numbers. These 
equations for R£m{t,r) simplify considerably in the adiabatic limit, in which 
Rem{t,r) oscillates much more rapidly than does the background spherically- 
symmetric geometry 7^^ . In this limit, the mode equation resembles the corre- 
sponding equation for a masslcss scalar field in an exact Schwarzschild geometry 
[10,11], except that the gravitational field and potential vary slowly with time. 
This allows one to approximate much of the classical space-time by means of 
the Vaidya metric [12]. Sec. 3 treats the boimdary conditions for suitable radial 
functions on the final surface T,p . It also begins the process of evaluating the 
second-linearised classical action S^l^^ corresponding to spin-0 perturbations 
on Sir , as needed in finding the quantum amplitude. In Sec. 4, the analytic be- 
haviour of S^^J^^ is studied in detail, as a function of complex T , leading to an 
expression (2.1,4.16) for the amplitude for purely spin-0 perturbations on Si? , 
found by taking the limit ^ — > 0+ , with T = \T\ exp{—i9). Sec. 5 is concerned, 
as a practical or efficiency measure, with the question of extracting standard 
dimensionful factors from all quantities appearing in the calculations, so that 
all equations can be understood more clearly, reduced to relations between di- 
mensionless quantities. Sec. 6 contains a short Conclusion. 



2 Adiabatic radial functions 

In this paper, we are mainly concerned with evaluating the scalar-field contribu- 
tion to quantum amplitudes for diffuse weak-field gravitational/scalar data on 
the final surface, which is separated from the initial surface by the 'Euclidean 
time-interval' r , or equivalently the 'Lorentzian time-interval' T , related by 
T = iT . Here, as in Eq. (1.2.3), we shall generally take T itself to be complex, 
of the form T = |r| exp(— i^?), where < 6* < tt/2, in order that the classical 
boundary-value problem should be well-posed. For such amplitudes, following 
Section 1.5, we need to compute 

Amplitude = (const.) x expji S'dass (/ijj, 0)/ ; (/iy, 0)f ; T | , (2.1) 

where hij = gtj denotes the intrinsic 3-metric on the bounding surface, and 
[Eq.(I.5.1)] 

Iff f\,, . 1 



(2.2) 

gives the Lorentzian classical action, apropriate to specifying hij and </> on both 
bounding surfaces, each surface having the same ADM mass M , where the 
Lorentzian time-separation between them is T . Equivalently, one may write 
Amplitude = (const.) x exp(— /dass), where /class = — ^'S'ciass is the 'Euclidean 
classical action'. 
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In the present context of ncarly-sphcrical collapse to a black hole, where 
the matter source is a massless scalar field </> , we expand the fields about the 
spherically-symmetric background solution (7^,^ , $), as: 

g^,u = 7/*. + hi'} + + . . . , (2.3) 
(j) = ^ + ^(1) + ^(2) + . . . . (2.4) 

The classical action then takes the perturbative form 

-Sclass = 'S'cliss + -^cliss + "^cliss + • • • ) (2-5) 

where S^f^ss background action, and the next-order correction S^^^^ is of 

second order in small perturbations. Explicitly [Eq.(I.5.10)], 

^ i (1 - /J ' \ (L-IJ '™ ■ 

(2.6) 

The quantity of physical interest is the amplitude corresponding to the weak- 
field non-spherical data {h^ij') p and ((/)(^))^ on the final surface, given here 

(approximately) by exp(iS'^[g|gg). For simplicity of exposition, we shall from now 
on assume that the initial data are chosen to be exactly spherically-symmetric, 
given only by (7^ , $)/ . Equivalently, we now take 

= 0, (0(1)), = 0. (2.7) 

Then, the amplitude cxp(zS'ciass) will only depend on the contributions at the 
final surface in Eq.(2.6), [which themselves depend on {h^^j'p , ; T)]. As a 

practical matter, we could easily put (/i^^^), back into the calculations that 
follow. Physically, the analogous step of 'turning back on the early-time pertur- 
bations' corresponds, in 'particle language' rather than in the 'field language' 
being used in this paper, to the inclusion of extra particles in the in-states, to- 
gether with the original spherical collapsing matter, and asking for the late-time 
consequences. This was first carried out by Wald [13]. 

In this paper, we concentrate on the scalar-field contribution to the quantum 
amplitude exp(iS'ciass)- That is, we need to compute 

^il, scalar = \ I 7r« <j>^'^ , (2.8) 

where the linearised perturbations (h^fl} ,(j)^'^) obey the linearised field equa- 
tions (1.3.20-22,25) about the spherically-symmetric background (7^1/ , <&). Here, 

{h^}} ,(j)^'^) must agree with the prescribed final data (hij^ (j)^'^) p at the fi- 
nal surface S p , and be zero at the initial surface E/ . In the Riemannian 

case, with a real Euclidean time-interval t between S/ and T,p , or in the case 
(1.2.3) of a complex time-interval T = t exp(— i0) between E/ and Yip , where 
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< < 7r/2, one expects that this hnear boundary- value problem will be 
well-posed. The other, gravitational, contribution 

^il„.av = 3^ /^^ .(^)- /.g) (2.9) 

will be studied subsequently. 

As described in Section 1.4, at late times the perturbed scalar field equation 
reduces to 

V'^V^0(i) = , (2.10) 

with respect to the spherically-symmetric background j^i, . Here, in contrast to 
the case of Sees. 1. 3. 4, it is most natural to describe the gravitational field in the 
Lorentzian case: 

ds^ = _ eb{t,r) de + e»(*''-) dr^ + {dO^ + sm^Od^^) . (2.11) 
Making the mode decomposition with respect to 'Lorentzian coordinates' {t, r, 6, (j)) : 

^ oo t 

cj>^^\t,r,e,4') = - ^ ^ Ye^{n) Ri^{t,r) , (2.12) 
one arrives at the {£ , m) mode equation: 

(e^b-"y^dryRem~{dtyRem-^(dt{a-b)){dtRim)~Viit,r)Rem = . 

(2.13) 

Here, 



where m{t ,r) is defined by 

2m{t, r) 



exp 



(-a(t,r)) = 1-^ . (2.15) 



For high frequencies of oscillation, as measured in the nearly- Lorentzian case, 
it becomes simpler to understand the solutions of the mode equation (2.13). 
Noting that the present {t , r)-coordinate description is obtained by replacing r 
by it , consider a solution Rem{t , r) of the form 

Remit, r) ~ exp{ikt) ^kem{t,r) , (2.16) 

where S,him{t, ^) varies 'slowly' with respect to t . In particular, at spatial infin- 
ity, where r — » oo , Remit, r) is required to reduce to a flat space-time separated 
solution in which £,kem{t,r) loses its t-dependcncc [see Eqs.(2.18,21) below]. 

We are studying the boundary-value problem for real scalar perturbations 
(j)^^^ , as functions of {t,r,9,(j)), or equivalently for real functions Remit, r) as 
in Eqs.(2.13-15), subject to the initial condition (j)^^^ \t=o = and a final 
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condition in which 0^^^ \t=T is prescribed, where again the final time T is of 
the form T = |T|exp(— for < < tt/2. For ease of visuahsation, one 
may regard 6 as being small and positive. Were the propagation simply in flat 
space-time, the solution would be of the form 

^^^^ = - E E * ^'^^rnir) YU^) , (2.17) 

where the {aMm} are real coefficients and each function ^kim{r) is proportional 
to a spherical Besscl function rje{kr) [14]. In our gravitational-collapse case, 
^kim becomes a function of t as well as of r, but otherwise the pattern remains: 



= ^ E E / 

1=0 m=-e •'~ 



sin(fci) 
sm(fcT) 



Here, the {akem} characterise the final data: they can be constructed from the 
given 0^^^ \t=T by inverting Eq.(2.18). The functions S,kemit,r) arc defined in 
the adiabatic or large-|A;| limit, as in the previous paragraph, via Eq.(2.16), 
where Remit, r) obeys the mode equation (2.13). 

More precisely, provided that k is large, in the sense that the adiabatic 
approximation 

|fc| » i|d-6| , (2.19) 

|fc| » ll^l , » ll^l , (2.20) 

holds, the mode equation reduces approximately to 

e(6-a)/2 d_ ^ _ ) ^^^^ ^ _ (2.21) 

Here, of course, the functions e^''"''^/^ and Ve do still vary with the time- 
coordinate t , but only adiabatically or 'slowly'. 

As described further in [2,12,17], the geometry in the space-time region is 
expected to be approximated very accurately by a Vaidya metric [12], corre- 
sponding to a luminosity in the radiated particles which varies slowly with time. 
Such a metric can be put in the diagonal form 

e-« = l-'JH^ ; e" = (jP^Xe- , (2.22) 

r \f{m)J 

where m(t,r) is a slowly-varying function, with m = (dm/dt), and where the 
function /(m) depends on the details of the radiation. Then Eq.(2.19) implies 
that 

777 

|fc| » l-l , (2.23) 

TO 
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provided that 2m{t,r) < r < 4m.{t,r) . In this case, the rate of ehange of the 
metric with time is slow compared to the typical frequencies of the radiation; 
further, the time- variation scale of the background space-time 'y^vix) is much 
greater than the period of the waves. With frequencies of magnitudes ^ 
dominating the radiation, and with |m| of order [15], the adiabatic 
approximation is equivalent to w? » 1 , which corresponds to the semi-classical 
approximation. If. as expected [15], is a measure of the time taken by the 
hole to evaporate, then r < Am <^ , provided that :s> 1 ■ Thus, in the 
large— A; approximation used in deriving Eq.(2.21), it is valid at lowest order to 
neglect time-derivatives of the background metric, out to radii small compared 
with the evaporation time of the hole and with the time since the hole was 
formed. 

It is natural to define a generalisation r* of the standard Regge- Wheeler 
coordinate rj for the Schwarzschild geometry [10,16], according to 

Under the above conditions, the time-dependence of r*{t, r) is negligibly small, 
and one has r* r* for large r , where 

r* = r + 2Mln((r/2M) - l) (2.25) 

is the Regge- Wheeler coordinate, expressed in terms of the Schwarzschild radial 
coordinate r . In terms of the variable r* , the approximate (adiabatic) mode 
equation (2.21) reads 



+ (fc^ - Vi) = . (2.26) 



3 Boundary conditions 

We consider here, for definiteness, a set of suitable radial functions {^kemir)} on 
the final surface Sj? . As above, since the mode equation (2.13) does not depend 
on the quantum number m, we may choose ^kem{r) = ^ke{f) , independently 
of TO. 

We seek a complete set, such that any smooth perturbation field (j)''^'^ (T, r, 9, cp) 
as in Eqs.(2.12,13), restricted to the final surface {t — T}, of rapid decay near 
spatial infinity, can be expanded in terms of the ^kemir)- The 'left' boundary 
condition on the radial functions {^keir)} is that of regularity at the origin 
{r = 0}: 

^€(0) = . (3.1) 
The solution to the radial equation which is regular near the origin is: 

^ke{r) = r^u{r) oc rj^ikr) oc | ((const.) x (fcr)^+i) -h 0((fcr)^+^) | , 

(3.2) 
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where, again, the are spherical Besscl functions [14], and we have assumed 
that, for small r , m(r) ~ , and have neglected O(r^) terms. These radial 
functions are purely real, for real k and r . For k purely real and positive, 
the radial functions describe standing waves, which, for mode time-dependence 
g±jA:<^ have equal amounts of 'ingoing' and 'outgoing' radiation. 

For the 'right' boundary condition, note that the potential Ve{r) vanishes 
sufficiently rapidly as r — > oo that a real solution to Eq.(2.26) behaves according 
to 

Cke{r) ~ (^Zke exp(ifcr*) + zl^ exp(-ifcr*)^ (3.3) 

as r — > 00 . Here the Zke are certain dimensionless complex coefficients, which 
can be determined via the differential equation by using the regularity at {r = 
0}. The (approximately) conserved Wronskian for Eq.(2.26), together with 
Eq.(3.3), and the property 

hm = iTvSik-k') , (3.4) 

r*~*co i^k — k') 

give the normalisation condition, for — oo < fc, /c' < oo and Roc, oo : 



r 

Jo 



dr e(»-'')/2 a^r) Ck'eir) = 2w {zuf (sik-k')+S{k+k')) . (3.5) 



This normalisation is only possible in an adiabatic approximation. Note that 
the radial functions {^fc^} form a complete set only for fc > , as a result of our 
boundary conditions. 

The above result makes it possible to evaluate the perturbative massless- 
scalar contribution to the total classical Lorentzian action ^dass = Sl^^l^^ + 
^dLs + ■■■ of Eq.(2.5), with S^f^ss given by Eq.(2.6). This contribution, namely 
^iiLs scalar Eq.(2.8), is then given in the notation of Eq.(2.12) by 



sSL[<t>^'^;T] = E / dre^'^-'y^Rem{dtR}j\^ , (3.6) 

e=0 m^-e •'O 



smce 

dVL Yim y^'m' = ^U' ^mm' ■ (3.7) 

Within the adiabatic approximation above, and using Eq.(3.5), this gives the 
frequency-space form of the classical action: 



o(2) 
class 



CO I -oo 

{akim};T = 7^ E E / dk k \zkt^ \akim+ a-kimt cot{kT) 

(3.8 

in terms of the final data {akem}- 
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From a mathematical point of view, one would expect to work only with the 
set of square-integrable scalar wave- functions on the final boundary S p , that 
is, the set L^{R^, dre^"-^'/^). To express this, define 

^em{r) = r J dn Ytmin) <p^'\t,r,fl)\ . (3.9) 

Then the square-integrability condition reads 

dr e'-'^-'^/^ \i,em{r)f < oo , (3.10) 

or, equivalently, 

/oo 
dk \zkef \akem + a-Mmf < oo . (3.11) 

The left-hand sides of Eqs.(3.10,ll) are in fact equal. This arises from the 
completeness property 



g(a-6)/2 r ikdr) ^kdr') ^ ^^sir-r') (3.12) 

-'-oo \Zki\ 

and the inverse of Eq.(2.18): 

aum + a-um = — -2 / dr e^"-^'^''^ ^,u{r) i^emir) ■ (3.13) 

2-K\zkA Jo 



From a physical point of view, one expects also that taking scalar boundary 
data which are not square-integrable will lead to various undesirable properties, 
such as infinite total energy of the system, or an infinite or ill-defined action. 



4 Analytic continuation 



The perturbative classical scalar action S^^^^^ of Eq.(3.8) was derived subject to 
the adiabatic approximation, and also to the requirement that the time-interval 
T between the initial and final surfaces, measured at spatial infinity, is complex, 
of the form T = \T\ exp{—i6), provided that < ^ < 7r/2 . In this case, the term 
A;cot(A;T) in the integrand of Eq.(3.8) remains bounded for < fc < oo , and 



;(2) 



{akem}', T 



given 



one expects to obtain a finite complex-valued action S^^^^^ 
square-integrable data 0'^' on the final surface Si? . Further, the dependence of 

(2) ' 

the complex function S^^^^ {(ikem}':T on the complex variable T is expected 

to be complex-analytic in this domain (0 < 6 < 7r/2), and, following Feynman 
[6], ordinary Lorentzian-signaturc quantum amplitudes should be given by the 
limiting behaviour of exp(iS'^^j]gg) as ^ ^ 0+ . 
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If, on the other hand, one restricts attention to the exactly Lorcntzian- 
signature case {9 = 0), then the integral in Eq.(3.8) will typically diverge, due 
to the simple poles on the real-frequency axis at 

k = K = ^ , (4.1) 

(n=l,2,... ). 

Clearly, in order to be able to treat this classical Einstein/massless-scalar 
boundary- value problem in a way which is analytically sensible, one has to allow 
T to be displaced into the complex as above, whether by a small or a large angle 9 
of rotation. The spherically-symmetric 'background' 4-geometry 7^1^ and scalar 
field <I> will typically be complex. One might hope that the general (asymmetric) 
boundary-value problem of this type would admit a unified mathematical de- 
scription, based on the (conjectured) property that the boundary- value problem 
is strongly eUiptic [7], up to gauge. 

Given that T = \T\ cxp(—iS) is slightly complex, with < (5 ^ 1 , consider 



an integral such as Eq.(3.8) for S^^^^ {akem};T 



Write this as 



where 



/•oo 

J = E / * f^rnik) COt(fcT) , (4.2) 

em 

fim{k) = irk \zm\^ \aMm + a-Mryi? ■ (4.3) 

There arc infinitely many simple poles of the integrand at k = fc„ [n = 
1,2,... ), just above the positive real fc-axis. We then deform the original con- 
tour C along the positive real fc-axis into three parts, Cf , and Cq, , where 
< a ^ 1 . The contour C, lies in the lower half-plane, half-encircling each 
of the simple poles near the positive real fc-axis, with radius e . The curve Cb. , 
also in the lower half- plane, is an arc of a circle |fc| = i? of large radius. The 
curve Ca is part of the radial line arg(fc) = —a . We write 



(4.4) 



Im +Cr-C^ 
= Ja + Jr + Je ■ 

Starting with the integral Jn , one finds 

\Jr\ < E / d9 R\fem{R,0)\ coth{\T\Rs-m9) , (4.5) 

em 

where k = i?e~*^ on Cr, and we have used |cot(A:T)| < coth(|T|i?sin0) . 
When the limit i? — > 00 is eventually taken, one expects that the contribution 
from Cr to the total action should vanish; this requires that \ftm{k)\ should 
decay at least as rapidly as |fc| , as |fc| ^ cx) . In fact, on dimensional grounds, 
one expects that 

\hmm ~ (4.6) 
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as — > 00. To see this, rewrite the radial equation (2.26) in terms of the 
operator 

£^ = ,(6-a)/2|.j^g(6-a)/2|.()^_^^(^) , (4.7) 

which is self-adjoint with respect to the inner product (3.5). Then note that 
Eq.(3.13) can be rewritten as 

aum + a-um = — — ^ / dr e^""'')/^ ^fc^(r) Cei^emir) ■ (4.8) 

ZTrfc^ \zkt.\ Jo 

We have used the boundary condition (3.1) and assumed that tpem{r) dies out at 
large r . The form (4.8) is just an expression of the self-adjointness of the radial 
equation. Now consider the dimensions of the quantities involved - treated 
explicitly in Sec. 5 below. In particular, iplmir) has dimensions of length and 
\zuf' is dimensionless. In the limit i?oo oo , and for large k (so taking a WKB 
approximation for the radial functions), the integral in Eq.(4.8) can only involve 
the dimensionless frequency 2Mk , where M is the total mass (true ADM mass) 
of the space-time. This gives the desired behaviour (4.6) at large |fc| . 

The contour gives a purely imaginary contribution to the total Lorontzian 
action; also (below) the curve Ca gives a complex contribution. We shall inter- 
pret the quantity cxp[— 2 Im(S')] , up to normalisation, as describing the con- 
ditional probability density over the final boundary data. To compute , we 
assume that fem{k) is analytic in a neighbourhood of k = dn , where 



= , (4.9) 



for n = 1,2,... . Then 



lim I dk fim{k) cot(fc|T|) 



(4.10) 



For the curve Ca , one has 



Jo. = - d\k\ e-"' hm(\k\,a) cot(\k\e-"'\T\) . (4.11) 

We shall need the properties [14] 

oo 

n=—oo ^ ' 

and 

i— — = p.p.^— ^ tT,5(x-a) , (4.13) 

[x — a±ie) [x — a) 
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where P.P. denotes the principal part. Assuming that fem{k) is regular along 



Ca , one has, for small a : 

JW = lim(l-ia)5: ± f\\k\ . (4.14) 

a^o+ ^ ^ Jo {\kT\ - n-K - la.) 

In the further limit i? — > oo , this gives 



= p-v. + ^ E E /^'»(^") ' (4.15) 

' ' em n=l 

where P.V. denotes the principal-value part of the integral. 

Using Eqs.(4.4,5,10,15), the classical action for masslcss scalar-field pertur- 
bations, in the case that T = \T\ ex.p{—i6) is very slightly complex, is 



c(2) 
class 



2 - oo t oc 

{aki7n};\T\ = realpart + E E E f^rni'^n) 

I I e=o m=-e n=l 



= realpart + E (^nlznif \anem + a-nimf ■ 

(4.16) 

The real part of S^^^ is, of course, also calculable from the equations above. It 

contains the principal- value term and the real part of Eq.(4.10). The main, semi- 
classical contribution to the quantum amplitude is then exp ^iS'^j^^^gg [{«fe£m}; |r|] ^ . 
The probability distribution for final configurations involves only Im(S'^^j]j,j,); 

(2) 

the more probable configurations will have S),^!^^ lying only infinitesimally 
in the upper half-plane. Whether probable or not, those final configurations 
{flfc^m} which contribute to the probability distribution must yield finite ex- 
pressions in the infinite sums over nl in Eq.(4.16). There will be a corre- 
sponding restriction when the data are instead described in terms of the spatial 
configurations {'0£m(r)}. Also, as can be seen in [17], the complex quantities 
Zni {p-ntra + Of-nim) appearing in Eq.(4.16) are related to Bogoliubov transforma- 
tions between initial and final states, thus providing a further characterisation 
of the finiteness of Im(5^j2ss) in Eq.(4.16). 

With regard to the sum over £ in Eq.(4.16), one imagines that a cut-off ^max 
can be provided by the radial equation (2.26). In the region where (V^(r) — fc^) > 
, one has exponentially growing radial functions, whereas for iyiir) — fc^) < , 
one has oscillatory radial functions. One defines £max by (V^max('') — fc^) = , 
and restricts attention mainly to oscillatory solutions. 

When one has both initial and final non-zero Dirichlet data labelled by 
'coordinates' {cife^^} and {a^^^}, the perturbative classical scalar action S'^^ss 
includes separate terms of the form (4.16) for the initial and final data. But 
•^cfass includes a cross-term between a'l}^ and a[,^^ , which represents the 
correlation or mixing between the initial and final data. The total action will 
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naturally be symmetric in aj^fl^ and ai.(„^ , and the coefficients z„£ will be the 
same (they are time-independent) up to a phase. For large \T\ , the cross-term 
becomes negligible, and one has two independent contributions to the classical 
action, one being a functional of {o^^^}, the other of {ctfef^}- 



5 Dimensional analysis 

In Sec. 4 above, we made use of a dimensional argument, in order to obtain the 
estimate (4.6) of the rate of fall-off of \fem{k)\ for large |A;| . In the present 
Section, we give a more thorough treatment of the dimensionality of the main 
quantities appearing in the preceding paper I, this Paper II and subsequent 
work. 

The classical perturbative Einstein equations are schematically of the form 
^ [dcj)^^^) , where L is a length-scale characteristic of the gravitational 
field. Hence, the massless-scalar perturbations i/)'^^^ are dimcnsionlcss. Similarly, 
one finds that the metric perturbations h^}} are dimensionless. From this, one 
deduces that the functions ^kem{tti~) of Eq.(2.16), appearing in the adiabatic 
approximation, are also dimensionless, as are the complex coefficients Zki of 
Eq.(3.3). 

In describing dimensionful quantities, let us denote by M the ADM mass 
of the 'space-time'. We recall again the distinction which it was necessary to 
draw in Sec. 1. 5 and in [3] between the naive mass, computed from the fall-off 
of the spatial metric gij of a hypersurface badly embedded near spatial infinity, 
and the true ADM mass. Here, M denotes the true ADM mass, as given by the 
standard definition for a hypersurface which approaches spatial infinity in the 
familiar way [10]. Suppose that the time-separation T at spatial infinity obeys 
\T\ ^ 2M , as will normally be the case; further, let the radius Rao tend to 
infinity. Then, in our treatment of massless perturbations, M can be regarded 
as setting a reference length-scale. 

We now write 

X = 2Mk (5.1) 

for the dimensionless frequency. Further, recall that the coordinates {akim} of 
Eq.(2.18) for the perturbative scalar field on the final surface T,p are defined 
irrespective of T. Prom Eq.(2.18), it follows that each akim scales as M^. That 
is, one can write 

akim = M'^ yem{x) , (5.2) 

for some (smooth) dimensionless complex functions yem{x)- This property can 
also be derived from Eq.(3.13). By following these arguments, one can write the 
imaginary part of the total classical action as 

^^{s2,sl{ytm};\T\]^ = (f) Yl (^^") 1^"^!' 12/"^™!' ' (5.3) 
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where 



TITT 

Xn = 2Man , " |T| ' '^^■^^ 

27rM , , 

= -jj^p • (5-5) 

Clearly, some understanding of the numerical magnitude of the s\im in 
Eq.(4.16) is necessary, if exp[— 2Im(5'^^j^gg)] is to be interpreted as a conditional 
probability density. (Indeed, the requirement that these probabilities sum to 1 is 
part of the requirement of unitarity!) In this lowest-order perturbative approx- 
imation, the probability density above for scalar fluctuations is Gaussian; the 
same holds in the present Einstein/massless-scalar case for the distribution of 
gravitational- wave fluctuations on the final surface E p ; this is treated in more 
detail in [18]. The Gaussian property is a perturbative aspect of a more general 
property - arising in the Positive Action conjecture [19,20]. For simplicity of ex- 
position, let us consider this for the case of Riemannian 4-metrics, although we 
expect some version to hold in the more general complex case. Suppose that, 
as in the 'Euclidean' path integral, we take all smooth 4-metrics g^j^^j , scalar 
fields (/),... which agree with Dirichlet boundary data /i^, (f>, . . . prescribed 
on a compact boundary dV, bounding a compact manifold-with-boundary V. 
The Riemannian action functional Ilg^v; ^; • • • ] may be seen to take arbitrarily 
negative values, when one applies suitable high-frequency conformal transfor- 
mations g^i, ^^g^iv to the metric [19]. But, at least for Einstein gravity with 
the simplest topology, the classical action /class is a non-negative functional of 
the boundary data [20] . Similar positivity properties can be investigated in the 
case of Einstein gravity coupled to matter. 

This dimensional analysis sets a natural normalisation for the (equal) quan- 
tities on the left-hand sides of Eqs.(3.10,ll). Using Eq.(5.2), one has 

/■oo /"OO 

2 ^ / dk \zkef \akem + a-kemf = ^ dx \ze{x)f \yem{x)f 

OC to ■ 

(5.6) 

The last line applies to semi-classical collapse and subsequent evaporation of 
a non-rotating black hole, where to oc is the time taken for complete 

evaporation [15]. 



6 Conclusion 

In this paper, we have derived the quantum amplitude [through Eq.(4.16)] for 
a spherically-symmetric configuration {hij,(p)i on the initial surface S/ to be- 
come a given configuration (/ly ,(j))p on the final surface S p , with Lorentzian 
time-interval T at spatial infinity, provided that the final 3-dimensional met- 
ric hijF is also spherically symmetric. In the amplitude, which is of the form 
(const.) X exp(iS'^^j^gg) , the action iS^f^ss depends approximately quadratically on 
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the (non-spherical) perturbative part of the final data (f)F ■ Further, S^^l^^ has 
both a real part and an imaginary part. The imaginary part leads to a Gaus- 
sian probability density oc exp(— 2 Im(S'^j^^^gg)) , while the real part gives 
rapid oscillations through the phase of the quantum amplitude or wave function 
Vt. Corresponding quantum amplitudes for spin-1 and spin-2 fields have been 
calculated [18], while the fermionic massless spin-i case is summarised in [21]. 

The familiar description of black-hole evaporation in terms of Bogoliubov 
coefficients occurs in many works [22-24]. Following the 'Complex Approach' 
and the 'Spin-0 Amplitude' for black-hole evaporation, we next make a con- 
nection between the description and calculation of quantum amplitudes in the 
present paper and the alternative description in terms of Bogoliubov coeffi- 
cients [17]. Further, we treat the approximation of the radiative part of the 
space-time by the Vaidya metric [12], related to the description of adiabatic 
scalar-field solutions in Sec. 2 of the present paper. Yet another description of 
quantum amplitudes will later be provided, relating Eq.(4.16) to coherent and 
squeezed states; this gives a more general conceptual framework with which we 
have already examined the amplitudes of the present paper and the higher-spin 
amplitudes [18,21]. 
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